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Comments  on  the  Colliding  Beam  Fusion  Reactor  Proposed  by  Rostoker, 
Binderbauer  and  Monkhorst  for  Use  with  the  p-nB  Fusion  Reaction 

1.  Introduction 


In  a  number  of  recent  publications1,2  and  progress  reports,3  Rostoker,  Binderbauer 
and  Monkhorst  (RBM)  have  advocated  a  colliding  beam  fusion  reactor  (CBFR)  for 
aneutronic  fusion  based  on  the  p-]1B  reaction.  This  scheme  was  originally  proposed  as  an 
energy  source  for  all-electric  naval  ships,  and  the  Office  of  Naval  Research  has  supported 
concept  studies  at  low  levels  for  two  years.  In  support  of  the  ONR  effort,  NRL  initiated  a 
small  effort  by  us  to  evaluate  the  p-nB  CBFR.  This  report  summarizes  our  conclusions, 
which  unfortunately  are  that  the  concept  is  fundamentally  flawed,  and  cannot  produce  net 
fusion  power. 

Since  the  beginning  of  the  fusion  program,  researchers  have  recognized  the 
tremendous  advantages  of  using  the  p-nB  fusion  reaction.  Since  all  of  the  fusion 
products  are  charged  a-particles,  it  is  widely  (but  not  universally4)  believed  that  a  p-nB 
reactor  would  be  free  of  the  radiation  problems  and  bulky  shielding  associated  with  the 
production  of  fast  neutrons,  and  fusion  output  energy  could  be  converted  directly  and 
efficiently  into  electricity.  The  fuel  would  be  cheap  and  abundant.  However,  there  are 
formidable  obstacles  to  the  use  of  this  reaction.  As  shown  in  Fig.  1,  the  reactivity  is  large 
only  if  the  relative  energy  between  a  colliding  proton  and  nB  nucleus  is  very  high,  in  the 
vicinity  of  Eo  =  600  keV,  where  there  is  a  resonant  peak  Ofo  in  the  fusion  cross-section 
Cp(E).  If  E  is  off  resonance  by  as  little  as  ±150  keV,  gf(E)  falls  to  about  half  of  oFo.  If 
the  protons  and  boron  are  thermally  distributed  at  a  temperature  of  650  keV,  the  reaction 
rate  is  only  about  a  third  of  what  it  is  if  every  p-nB  encounter  occurs  at  exactly  energy 
E0.  Since  the  fuel  must  be  so  energetic,  and  the  energy  output  is  only  8.7  MeV  per  fusion 
reaction,  the  output  fusion  power  (before  conversion  to  electricity)  can  be  no  more  than 
about  15  times  the  power  invested  in  injected  protons,  even  if  proton  acceleration  is 
100%  efficient  and  every  injected  proton  undergoes  fusion.  In  addition,  bremsstrahlung 
energy  losses  are  a  severe  problem,  because  of  the  presence  of  nuclei  with  Z=5.  When  all 
of  these  factors,  and  others,  are  taken  into  account,  the  generally  accepted  conclusion5  (as 
restated  by  RBM1)  is  that: 

“For  a  thermal  p-^B  reactor,  the  electromagnetic  radiation  energy  is  greater  than  the 
nuclear  energy  produced,  and  a  reactor  that  produces  net  energy  is  possible  only  if 
the  conversion  efficiency  is  nearly  100%.” 

In  order  to  improve  the  energetics  to  the  point  where  net  energy  gain  might  be 
possible,  the  primary  objective  of  the  RBM  plan  is  to  have  nearly  every  encounter 


Manuscript  approved  September  30, 1998. 


1 


between  a  proton  and  a  nB  nucleus  occur  at  relative  energy  close  to  Eo.1"3'6"7  To  bring 
this  about,  they  plan  to  inject,  confine,  and  maintain  the  protons  and  nB  as  separate  cold 
beams:6 

“In  order  to  exploit  resonance  in  a  steady-state  reactor,  it  is  necessary  to 
maintain  the  proton  and  boron  beams  at  an  average  energy  difference  of  580  keV. 

In  addition,  the  temperatures  of  the  beams  must  be  substantially  less  than  140  keV.” 

Specifically,  RBM  propose  to  confine  the  reactants  within  a  field-reversed  configuration, 
as  shown  in  Fig.  2,  with  the  protons  maintained  as  a  rapidly  rotating  cold  beam  with  a 
thin  annular  spatial  profile.  The  nB  ions  form  a  fully  stripped  cold  species,  which  may 
be  a  stationary  component  or  a  slowly  rotating  beam.  The  relative  energy  of  the  p  and  nB 
is  Eo,  but  the  temperatures  of  each  species  must  be  no  more  than  tens  of  keV.6  The 
electrons  form  a  third  species,  which  rotates  at  a  speed  intermediate  between  the  p  and 
nB  beams.  Because  of  the  electrical  current  which  is  carried  primarily  by  the  rotating 
protons,  the  magnetic  field  reverses  in  crossing  the  ion  ring.  In  various  design  options, 
RBM  choose  the  magnetic  swing  to  be  in  the  regime  of  100  to  200  kG,  which  implies  that 
the  product  of  the  proton  density  np  and  the  thickness  Ar  of  the  ring  is  npAr  =  fewx  1015 

-J? 

cm  .  The  radius  of  the  proton  ring  is  taken  to  be  anywhere  from  30  cm  to  80  cm,  and  the 
thickness  of  the  ring  is  assumed  to  be  a  few  cm.  Since  the  equilibrium  is  expected  to  be 
at  least  approximately  of  rigid  rotor  type,  i.e.  the  ion  rotational  velocity  is  proportional  to 
the  radius  r,  it  is  important  that  the  proton  ring  be  thin  to  exploit  the  resonance  in  Gf . 

Forty  years’  experience  provides  convincing  proof  that  fusion  is  a  much  more 
challenging  problem  than  was  envisioned  by  the  pioneers.  Many  plasma  confinement 
schemes  have  been  studied.  In  considering  a  new  confinement  scheme,  the  standard 
procedure  is  to  consider  first  the  “classical”  confinement  and  transport  properties,  i.e. 
those  resulting  from  ordinary  Coulomb  collisions  between  the  charged  particles.  These 
are  relatively  straightforward  to  understand,  although  detailed  calculations  may  be 
difficult  in  complex  confinement  geometries.  Normally,  the  only  confinement  schemes 
which  are  considered  for  further  theoretical  and  experimental  study  are  those  for  which 
the  classical  properties  are  found  to  be  satisfactory.  Typically,  these  more  detailed 
investigations  uncover  a  variety  of  instabilities  and  collective  effects  which  may  interfere 
with  the  achievement  of  adequate  particle  and/or  energy  confinement,  and  which  must  be 
overcome.  However,  the  p-nB  CBFR  proposed  by  RBM  fails  to  pass  the  first  test. 
Simple  classical  collisional  properties  and  basic  physical  principles  indicate  that  the 
equilibrium  required  for  the  CBFR  cannot  be  sustained  for  long  enough  to  provide  fusion 
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gain.  Various  aspects  of  these  problems  have  been  elucidated  in  comments  sent  to 
Science  magazine  by  seven  scientists  in  response  to  Ref.  1.  "  All  of  these  authors  note 
that  the  p-nB  fusion  reaction  is  at  best  marginal  as  regards  energy  gain,  and  that 
everything  must  break  just  right  if  the  process  is  to  have  any  chance.  References  10-13 
also  present  a  variety  of  collisionality  arguments  similar  to  those  presented  here. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  2,  we  summarize  our  notation. 
In  Sec.  3,  we  discuss  in  some  detail  a  variety  of  collision  processes  which  destroy  the 
CBFR  equilibrium,  on  time  scales  much  faster  than  fusion.  In  Sec.  4,  we  present  an 
explicit  calculation  of  the  effect  of  p-nB  collisions  on  the  proton  temperature.  This  has 
been  an  ongoing  point  of  contention  between  RBM  and  ourselves,  and  we  also  point  out 
an  error  in  the  RBM  analysis  which  led  them  to  conclude  that  the  protons  can  be  cooled 
by  p-uB  collisions.  In  Sec.  5,  we  briefly  discuss  the  state  of  the  art  in  beam  source 
technology;  indications  are  that  the  beam  source  for  a  p-nB  CBFR  reactor  would  be  too 
large  to  meet  the  Navy’s  need  for  shipboard  power.  We  also  offer  brief  comments  on 
instability  concerns,  and  on  the  limited  experimental  data  base  for  field-reversing  ion 
rings.  In  Sec.  6  we  summarize  our  conclusions.  Detailed  Fokker-Planck  calculations, 
which  are  the  basis  for  Secs.  3  and  4,  are  deferred  to  the  Appendix. 

In  drafting  this  paper,  our  intention  has  been  to  make  the  arguments  accessible  to 
the  physics  community  at  large,  and  not  only  to  those  who  have  a  detailed  background  in 
plasma  physics. 

2.  Summary  of  Notations 

To  facilitate  comparisons  between  our  work  and  the  published  papers  of  RBM,  we 
have  tried  to  use  similar  notation;  however,  in  some  cases  we  have  chosen  notations  that 
are  a  bit  clearer  or  more  intuitive  to  us.  The  definition  of  quantities  such  as  Vj,  Vj  and  Tj, 
is  identical  to  the  usage  of  RBM.  We  use  cgs  units. 

Species:  Subscripts  p,  B  and  e  represent  protons,  nB  nuclei  and  electrons,  respectively. 

Macroscopic  parameters  and  distribution  functions: 

•  nj  is  the  number  density  of  species  j  at  a  particular  spatial  location;  j  =  p,  B,  or  e. 

•  pj  =  mjnj  is  the  mass  density  of  species  j  at  a  particular  spatial  location. 

•  Vj  is  the  mean  or  fluid  velocity  of  species  j  at  a  particular  spatial  location,  defined  as 

Vj  =  Jd3vvfj(v).  (1) 
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•  fj(v)  is  the  velocity  distribution  function  of  species  j,  normalized  to  unity. 

•  Vj  is  the  magnitude  of  the  vector  velocity  Vj. 

•  Tj  is  the  temperature  of  species  j,  defined  as 

Tj=jd3vfmj(v-Vj)2fj(v)-  <2> 

The  definition  of  temperature  does  not  assume  that  the  distribution  fj(v)  is  Maxwellian, 
or  even  that  it  is  isotropic,  but  in  thermal  equilibrium  it  reduces  to  the  usual  definition. 

.  Vj  =  ^Tj/nij  is  the  thermal  velocity  of  species  j. 

•  Ar  is  the  width  of  the  annular  proton  beam. 


Cross-sections,  rates,  and  characteristic  times: 

.  gf(E)  is  the  p-nB  fusion  cross-section,  as  a  function  of  E,  the  relative  energy  between  a 
colliding  proton  and  nB  nucleus.  Gfo^f<Eo)  is  the  peak  value  of  Of(E).  RBM15  cite 
two  different  references  for  Gf(E).  Feldbacher  and  Heindler  (1988)16  give  Gfo  = 
0.8xl0~24  at  Eo  =  620  keV,  while  Becker,  Rolfs  and  Trautvetter  (1987)17  give  Gfo  = 
1.15xl0-24  at  Eo  =  580  keV.  In  addition,  RBM  argue18  that  it  may  be  possible  to 
increase  Gf(E)  by  a  factor  of  1.6  by  spin-polarizing  the  fuel  nuclei.  Throughout  this 
paper,  we  shall  use  the  most  optimistic  estimate  of  Gf(E),  i.e.  the  Becker  et  al  value, 
with  the  assumption  of  100%  spin  polarization.  Thus,  we  assume  Gfo  =  1.8x10  24,  and 
we  round  off  Eo  to  -600  keV.  None  of  our  conclusions  are  sensitive  to  the  exact  value 
of  Gf(E)  or  Eo- 

.  Vf = nBv0GF(Eo)  is  the  peak  rate  at  which  a  proton  undergoes  fusion.  Here,  v0  =  1 .  lxlO9 
cm/sec  is  the  proton  velocity  at  energy  Eo-  This  rate  also  assumes  fuel  that  is  optimally 
spin-polarized  for  fusion. 

•  Vjk  is  the  effective  momentum  exchange  collision  frequency  for  species  j,  due  to 
multiple  Coulomb  collisions  with  particles  of  species  k. 

•  Vjk,E  is  the  rate  of  energy  change  for  species  j  due  to  Coulomb-collisional  energy 
exchange  with  species  k. 

•  t= v-1  is  the  characteristic  time  for  any  process. 

3.  Collisional  Processes 


RBM  have  commented  that  a  complete  self-consistent  analysis  of  the  collisional 
dynamics  of  the  CBFR  would  be  a  difficult  computational  task,  given  all  of  the 
collisional  processes,  driving  forces,  sources,  sinks,  and  the  complexities  of  the  magnetic 
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confinement  geometry.  We  agree  with  this  assessment.  However,  it  is  relatively  simple 
to  perform  a  rigorous  analysis  of  each  of  the  collisional  effects  separately.  We  present 
just  such  an  analysis  in  this  section.  Since  RBM  assume  that  the  equilibrium  can  be 
maintained  at  least  long  enough  to  allow  fusion  of  the  protons  resident  in  the  system  (in 
fact,  their  ultimate  hope  is  to  maintain  a  steady  state),19  we  shall  assume  that  in  each 
scattering  process,  the  scatterers  present  a  steady  state  distribution.  In  particular,  we  shall 
begin  with  RBM’s  assumption6  that  the  proton  and  nB  constituents  each  form  a  cool 
beam,  with  relative  energy  -600  keV  between  the  two  beams.  We  shall  conclude  that  an 
equilibrium  of  the  type  proposed  by  RBM  cannot  be  sustained  for  a  time  long  enough  to 
allow  fusion  gain,  because  any  one  of  several  collisional  processes  dissipate  the  assumed 
CBFR  equilibrium  on  much  shorter  time  scales. 

A.  Time  Scales  for  Fusion  and  Collisional  Processes 


Table  1.  Proton  collision  rates 


Symbol 

Description 

Rate  (sec'1) 

Vf 

Fusion  rate  at  resonance  (polarized  fuel) 

2.0xlCT15nB 

VpB 

p-B  momentum  exchange  collisions 

7.4x1 0_14nB 

vpB,E 

p-B  energy  exchange  collisions 

1.4xlO-14nB 

Vpp 

p-p  collisions  (momentum  or  energy  exchange) 

7.2xl0_13np  (10  keV/Tp) 3/2 

3* _ 

Stopping  of  protons  by  p-e  collisions 

2.6x1 0_14ne  (10  keV/Te)3/2 

Table  2.  Boron  collision  rates 

Symbol 

Description 

Rate  (sec'1) 

vBp 

Acceleration  of  boron  by  B-p  collisions 

6.7xl0_15np 

vBe 

Acceleration  of  boron  by  B-e  collisions 

5.9xl0“14ne(10  keWTe)3/2 

The  rates  for  protons  to  undergo  various  collisional  processes20  are  shown  in 
Table  1.  Note  that  densities  are  in  units  of  cm-3.  We  notice  immediately  that  the  fusion 
reaction  rate  Vf  is  much  slower  than  the  relevant  rates  for  dissipative  collisional 
processes.  This  is  a  red  flag  -  especially  for  a  fusion  scheme  that  depends  on  colliding 
beams.  For  the  p-nB  reaction,  the  peak  fusion  rate  (per  proton,  using  the  most  optimistic 
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estimate  of  the  cross-section,17  and  assuming  100%  polarized  fuel  at  precisely  the 
optimum  collision  energy  Eo  =  600  keV  for  fusion)  is  Vf  =  nBctFoVo  =  2.0x10  15nB  sec  \ 
where  Ofo  is  the  peak  fusion  cross-section  and  vo=l.lxl09  cm/sec  is  the  relative  p-nB 
velocity  at  energy  Eo.  However,  the  collision  frequency  for  p-uB  momentum  exchange 
scattering,  i.e.  for  90°  scattering  of  the  protons,  is20  vpb  =  7.4x10  14nB  sec  *,  37  times 
faster  than  the  fusion  rate  Vf-  Momentum-exchange  scattering  has  the  primary  effect  of 
isotropizing  the  proton  beam.  This  destroys  its  assumed  beamlike  velocity  distribution, 
and  leads  to  spatial  spreading  so  that  the  protons  are  no  longer  confined  to  a  tight  annulus 
as  assumed  by  RBM.  Momentum-exchange  p-nB  collisions  also  exert  a  drag  on  the 
boron  beam,  with  the  effective  collision  frequency  Vbp  =  (npmp/nBmB)vPB  shown  in  Table 
2.  The  consequent  acceleration  of  the  boron  beam,  which  is  over  an  order  of  magnitude 
(depending  on  the  value  of  np/nB)  faster  than  fusion,  detunes  the  assumed  resonance  at 
energy  Eo.  Proton  energy  loss  due  to  p-nB  collisions  proceeds  at  a  ratevPB,E  = 
1.4xlO'14nB  sec*1,  slower  than  vpB  by  about  twice  the  mass  ratio  mp/mB,  but  still  faster 
than  Vf  by  a  factor  ~  7.  This  process  slows  down  the  protons,  thereby  additionally 
detuning  the  assumed  resonance  at  energy  Eo.  The  effects  of  p-nB  Coulomb  scattering 
on  the  protons  scale  in  the  same  way  as  the  fusion  rate  (proportional  to  nB),  and  thus  these 
rates  are  faster  than  the  fusion  rate  by  the  same  large  factor,  irrespective  of  the  choice  of 
the  density  ratio  np/nB.  In  addition  to  the  p-nB  scattering  processes,  there  are  various 
dissipative  processes  involving  collisions  between  other  species.  The  proton-proton 
Coulomb  scattering  frequency  is20  vpp  =  7.2xl0"13np  (10  keV/Tp)3/2.  For  all  scenarios 
considered  by  RBM  for  CBFR,  vpp  »  vF.  Proton-proton  collisions  thermalize  the  protons 
in  their  own  frame  of  reference,  i.e.  these  collisions  convert  a  spread  in  proton  momenta 
(which  results  from  p-nB  momentum  scattering)  into  a  proton  energy  spread,  further 
degrading  the  fusion  resonance.  Proton-electron  and  boron-electron  collisions  are  also 
faster  than  the  fusion  rate,  for  typical  CBFR  scenarios.  These  collisions  indirectly  couple 
the  protons  to  the  boron,  leading  to  a  friction  that  can  stop  the  relative  streaming  of  the 
ion  beams,  heat  the  electrons,  and  also  influence  the  ion  beam  temperatures. 

We  shall  now  discuss  these  collisional  processes  in  more  detail.  We  also  discuss 
the  possibility  of  introducing  driving  forces,  sources  and  sinks  that  can  overcome  the 
dissipation  and  sustain  the  equilibrium. 

B.  Proton-nB  Coulomb  Scattering 

The  velocity  distribution  of  p  and  nB  is  illustrated  in  Fig.  3a,  for  the  typical 
“colliding-beam”  equilibrium  envisioned  for  the  CBFR.6  In  Fig.  3,  the  horizontal  axis 
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schematically  represents  the  azimuthal  velocity  v<p,  which  is  the  direction  of  beam  flow. 
The  vertical  axis  in  this  2-D  picture  represents  the  two  dimensions  normal  to  the  flow,  vr 
and  v2.  The  protons  and  the  nB  are  each  assumed  to  constitute  a  cold  beam,  i.e.  the 
spread  in  velocities  of  each  component  is  small  compared  to  the  relative  azimuthal 
velocity  between  the  two  beams. 

In  this  subsection,  we  consider  the  effects  of  p-nB  scattering  on  the  proton 
distribution,  and  to  do  this  we  shall  construct  our  picture  in  the  frame  of  reference  where 
the  boron  component  is  stationary,  i.e.  we  shall  use  the  variable  u  =  v-Vb  to  specify  a 
proton’s  velocity.  To  simplify  the  discussion,  let  us  assume  that  the  boron  beam  actually 
has  zero  temperature.  This  is  the  most  favorable  situation  for  the  CBFR,  and  finite  boron 
temperature  only  makes  matters  worse.  The  nature  of  the  scattering  process  is  easy  to 
understand  intuitively.  Because  the  mass  ratio  mB/mp=ll  is  large,  the  primary  effect  on 
the  protons  of  multiple  Coulomb  collisions  with  the  nB  is  to  scatter  the  protons  through  a 
substantial  angle.  Energy  transfer  from  the  protons  to  the  nB  is  slower  by  a  factor 
2mp/mB  =  0.18.  Thus,  to  lowest  order  in  the  mass  ratio,  the  Fokker-Planck  equation 
reduces  to  simply  diffusion  of  the  proton  velocity  on  a  constant-u  spherical  shell  in 
velocity  space,  as  shown  in  Fig.  3b.  To  the  next  order  in  2mp/mB,  the  energy  transfer 
between  species  results  in  a  friction  that  draws  the  velocity  of  each  proton  toward  the 
mean  boron  velocity  Vb. 

We  prove  in  the  Appendix  that  this  picture  is  rigorously  correct:  for  scattering  of 
protons  off  any  cold  species  k,  the  Fokker-Planck  equation  reduces  exactly  to  diffusion 

on  a  constant-energy  sphere,  plus  friction;  there  is  no  diffusion  in  the  proton  energy 
1  2 

jmpu  .  According  to  Eqs.  (A10,14),  the  Fokker-Planck  equation  for  p-k  scattering  can 
be  written  in  spherical  velocity  coordinates  (u,0,\]/)  as 
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where  the  dynamical  friction  in  spherical  coordinates  is 
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the  transverse  diffusion  coefficient  is 
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and  X.  is  the  Coulomb  logarithm.  Note  that  these  velocity  coordinates  are  chosen  so  that 
the  origin  u  =  0  is  at  the  scatterer  (i.e.  boron)  velocity,  and  the  polar  axis  is  along  the 
difference  in  the  beam  fluid  velocities,  Vp-V k> 

The  relative  magnitude  of  the  friction  and  transverse  diffusion  coefficients 
depends  on  the  mass  ratio  mp/mk.  For  scattering  of  a  light  species  (e.g.  protons)  off  a 
heavy  species  (e.g.  UB),  Eqs.  (4,5)  show  that  transverse  diffusion  dominates,  as  stated 
above.  If  (for  the  moment)  we  neglect  the  friction  term,  the  solution  to  Eq.  (3)  may  be 
written  in  the  form 


f(v,t)  = 


S(u-Uq) 

47tv2 


g(0). 


(6) 


where  the  cold  beam  initial  condition  is  g(0,O)  =  28(02)  and  g(0,t)  is  the  solution  to 


3g  D ,  1  3  .  3 

— ^  - - sin0  — g 

3t  2u2  sin 0  30  30 6 


(7) 


Thus  the  proton  distribution  diffuses  on  the  surface  of  the  constant-u  sphere,  over  a  time 
scale  xpb  =  vPB_1,  as  shown  in  Fig.  3b.  Note  that  this  diffusion  process  steadily  increases 
the  proton  temperature  Tp,  as  defined  in  Eq.  (2),  but  nonetheless  all  protons  have  the  same 

energy  ^mpu2,  as  seen  in  the  boron  frame  of  reference.  In  an  anisotropic,  non- 

Maxwellian  plasma,  characterizing  a  species  by  a  single  “temperature”  does  not 
necessarily  convey  the  same  meaning  as  is  usual  in  thermodynamics. 

Thus,  if  p-nB  collisions  were  the  only  consideration,  the  constant-u  spherical 
shell  would  fill  in  completely,  as  shown  in  Fig.  3c,  at  a  time  of  order  2tpb  which  is  less 
than  3%  of  the  time  Tf=1/vf  required  for  fusion.  The  resulting  isotropic  distribution 
function  would  no  longer  be  a  “beam”  distribution,  and  it  is  incompatible  with  a  thin- 
annulus  proton  equilibrium,  or  indeed  with  any  rapidly  rotating  rigid-rotor-type 
distribution.  Thus,  the  p-nB  scattering  process  is  lethal  to  an  equilibrium  of  the  type 
assumed  by  RBM.  However,  the  situation  is  made  even  worse  by  a  peculiarity  of  the 
CBFR  field-reversed  configuration.  In  this  magnetic  confinement  geometry,  protons  are 
confined  only  if  they  rotate  azimuthally  in  the  “diamagnetic  direction.”  The  coordinates 
may  be  defined  so  that  this  is  the  direction  with  v<p>  0.  Protons  rotating  in  the  opposite 
direction  (v9  <  0)  promptly  escape  from  the  ends  of  the  device.21  As  the  protons  diffuse  to 
fill  the  energy  shell,  any  protons  scattered  into  the  left  hemisphere  in  Figs.  3b,c  (i.e 
negative  v9  or  0  >  %)  are  lost.  Thus  the  diffusion  process  does  not  actually  proceed  to  the 
point  of  uniformly  filling  the  energy  shell.  Instead,  about  97%  of  the  protons  are  lost  by 
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scattering  into  the  negative-v^  hemisphere  before  fusion  takes  place.  Since  the  energy 
invested  in  proton  acceleration  (taking  into  account  the  imperfect  efficiency  of  proton 
acceleration)  amounts  to  at  least  10%  of  the  fusion  output  energy  if  every  proton  fuses, 
this  in  itself  makes  it  clear  that  the  CBFR  cannot  provide  a  net  fusion  gain. 

C.  Effect  of  nB-p  Collisions  on  the  Boron  Beam 

We  have  previously  offered  the  argument  of  Sec.  3B  in  a  comment  on  Ref.  1 
which  was  sent  to  Science. 13  RBM  responded  as  follows:22 

“These  conclusions  are  inconsistent  with  conservation  of  momentum  and  energy 
that  must  be  satisfied.  Due  to  collisions,  protons  slow  down  and  boron  speeds  up. 
From  the  conservation  of  momentum 

npmpVp  +  nBmBVB=  (npmp+  nBmB)Vc, 

with 

- rynp -  - Mk - Vb  =  Vp 

npmp  +  nBmB  p  npmp  +  nBmB  p 

because  nBmB  /  npmp «  1  and  VB  «  Vp.  The  final  velocity  of  protons  is  only  a  little 
less  than  the  initial  value.  .  .The  protons  ‘isotropize’  as  noted  by  Lampe  and 
Manheimer,  but  in  a  frame  of  reference  moving  with  the  large  azimuthal  velocity  Vc 
so  that  the  number  of  protons  lost  because  they  rotate  in  the  wrong  direction  is 
negligible.  The  energy  loss  of  protons  is  much  less  than  600  keV.” 

In  essence,  RBM  object  to  our  treatment  of  the  boron  as  a  steady-state  scatterer 
species  with  a  constant  fluid  velocity  VB.  To  clarify  matters,  we  restate  the  thinking 
underlying  our  calculations.  The  calculations  of  RBM  all  assume  -  indeed  the  CBFR 
concept  only  makes  sense  for  p-nB  if  one  assumes  -  that  there  is  an  equilibrium  in  which 
the  p  and  the  nB  each  form  a  cold  beam,  with  relative  energy  ~  600  keV  between  the 
beams.  Our  objective  is  to  test  whether  this  assumed  steady  state  can  persist  for  a  fusion 
time.  Since  it  is  not  easy  to  do  a  complete  self-consistent  calculation  of  the  evolution  of 
all  of  the  species,  including  all  of  the  scattering  processes,  we  begin  by  accepting  the 
steady  state  assumption,  as  regards  the  scatterers  in  each  scattering  process,  and  address 
the  question:  If  the  scatterers  constitute  the  required  cold  beam,  is  it  possible  for  the 
scattered  species  to  maintain  its  cold-beam  distribution?  In  general,  the  answer  is  that  it 
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is  not.  The  required  steady  state  is  destroyed  by  whichever  process  turns  out  to  be  the 
fastest. 

With  this  preamble,  let  us  consider  the  effect  of  nB-p  collisions  on  the  boron. 
Our  Eqs.  (A10)  show  that  the  effect  of  p-nB  collisions  on  the  boron  is  primarily  a 
frictional  drag;  diffusion  (i.e.  heating)  of  the  nB  is  slower  by  the  mass  ratio  mp/mB.  The 
time  scale  for  friction  to  accelerate  the  boron  to  the  proton  beam  velocity  is  (by 
conservation  of  momentum)  xBp  =  tpB  pB/pp,  where  pB  =  nBmB  is  the  boron  mass  density 
and  pp  =  npmp  is  the  proton  mass  density.  In  much  of  the  work  of  RBM,  the  parameters 
are  chosen  so  that  pB  >pp,  in  which  case  proton  scattering  is  faster  than  boron 
acceleration.  But  for  the  parameters  cited  cited  in  Ref.  1,  np=400nB  and  therefore  pp  » 
pB.  In  this  situation,  acceleration  of  the  boron  beam  is  the  faster  process,  as  stated  by 
RBM.  The  center-of-mass  velocity  Vc  is  indeed  close  to  the  initial  proton  velocity  Vp, 
and  the  system  will  evolve  to  a  state  where  VB  =  VP  =  VC  over  the  time  scale  xBp  = 
1.5xl014np_1  =  3.7x10' ‘ns-1,  which  is  over  three  orders  of  magnitude  faster  than  the 
proton  fusion  time  Tf.  Thus,  for  this  particular  choice  of  parameters,  we  agree  with  the 
conclusion  stated  by  RBM  that  the  protons  will  not  end  up  with  a  large  velocity  spread. 
However,  the  consequences  of  boron  acceleration  are  even  more  devastating  for  the 
CBFR  concept:  the  protons  and  nB  end  up  as  co-rotating  beams,  rather  than  colliding 
beams.  The  temperature  of  each  of  the  beams  is  small  compared  to  the  600  keV  required 
for  efficient  p-nB  fusion  (see  Fig.  1),  and  virtually  all  encounters  between  protons  and 
nB  nuclei  occur  with  relative  energy  «  600  keV.  Consequently,  there  will  be  virtually 
no  fusion  reactions  at  all.  For  the  CBFR  to  make  any  sense  at  all6  for  p-nB,  there  must 
be  some  mechanism  (as  yet  unspecified)  that  maintains  the  protons  and  boron  as  separate 
colliding  beams  with  relative  energy  close  to  600  keV.  If  the  boron  is  somehow 
maintained  at  a  constant  rotation  velocity,  then  our  previous  argument  holds,  and  proton 
diffusion  leads  to  isotropization  about  the  boron  velocity. 

D.  Energy  Loss  Effects  Due  to  p-nB  Scattering 

In  the  previous  section,  we  have  seen  that  the  frictional  acceleration  of  the  boron 
beam,  due  to  p-nB  collisions,  can  be  important  when  pp  >  pB.  Up  to  this  point,  we  have 
neglected  the  energy  loss  of  the  protons  due  to  p-nB  scattering,  as  it  is  less  rapid  than  the 
pitch-angle  scattering.  However,  if  it  is  assumed  that  VB  is  maintained  at  a  constant 
value,  the  exact  solution  of  the  Fokker-Planck  equation  (3)  for  the  proton  velocity 
distribution,  including  both  the  friction  and  diffusion  terms  due  to  p-nB  scattering,  is  of 
the  form 
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(8) 


fj(u,e,q>,t) 


1 

4rcv2 


8(u-u0(t))g(6,t), 


where  uo(t)  is  given  by 


dug 

dt 


(9) 


F  is  the  friction  coefficient  of  Eq.  (4),  and  g(0,t)  is  the  solution  to  Eq.  (A18)  for  diffusion 
on  the  collapsing  spherical  shell.  Because  F  is  small,  g  is  nearly  identical  to  the 
diffusion  Green’s  function  defined  in  Eqs.  (6,7).  Thus,  the  proton  distribution  spreads  out 
diffusively  on  the  spherical  shell,  while  the  proton  energy  slowly  decreases  due  to 
friction,  i.e.  the  shell  as  a  whole  collapses  down  self-similarly  toward  Vb  because  of  the 
friction.  At  any  given  time  all  protons  have  the  same  energy,  i.e.  the  shell  does  not 
thicken,  as  there  is  no  diffusion  in  proton  energy.  (Recall  that  this  is  a  consequence  of  the 
assumption  that  the  boron  is  cold).  The  time  scale  for  the  loss  of  all  of  the  proton  energy 
is  vpb,e-1  ~  (mB/2mp)vPB_1,  which  is  still  about  seven  times  faster  than  the  fusion 
timescale.  Since  the  energy  loss  rate  is  rapid  compared  to  the  fusion  rate,  the  relative 
energy  between  the  protons  and  boron  cannot  be  maintained  at  the  resonant  value  600 
keV  which  is  necessary  for  efficient  p-nB  fusion.  Indeed,  the  fusion  rate  falls  off 
substantially  when  the  protons  lose  only  a  quarter  of  their  energy. 

Once  the  angular  direction  of  the  proton  velocities  has  been  spread  out  on  the 
energy  shell,  it  is  difficult  to  imagine  any  mechanism  that  could  reaccelerate  the  protons 
to  maintain  the  resonant  energy.  Any  force  that  acts  equally  on  all  protons  (e.g.  an 
electric  field),  will  merely  displace  the  spherical  shell  as  a  whole,  rather  than  expanding 
the  shell  as  would  be  necessary  to  restore  the  resonance. 


E.  Maintenance  of  the  Velocity  Separation  Between  p  and  nB  Beams 

Irrespective  of  whether  pB/pP  is  large  or  small,  p-nB  friction  brings  Vp  and  Vb 
together  in  a  time  small  compared  to  tf.  If  there  are  to  be  “colliding  beams,”  there  must 
be  other  forces  that  maintain  the  average  velocities  of  each  of  the  two  beams.  It  is 
difficult  to  see  how  this  could  be  done,  but  let  us  assume  there  is  some  such  force  active 
on  the  protons.  The  situation  is  then  analogous  to  ordinary  Ohmic  heating,  in  which  an 
electric  field  acts  on  electrons  to  maintain  a  constant  electric  current,  despite  friction  due 
to  electron-proton  pitch-angle  scattering.  In  the  present  case,  the  protons  are  the  light 
species,  and  boron  is  the  heavy  species.  As  in  ordinary  Ohmic  heating,  the  force  must 
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supply  energy  to  the  protons  at  a  rate  ^v^npHi^Vp^  -  VB,)  ,  which  is  about  four  times 

greater  than  the  fusion  power.  This  energy  is  dissipated  as  proton  temperature,  i.e. 
through  the  transverse  velocity  diffusion  process  discussed  in  the  previous  section.  So 
for  this  reason  as  well,  the  CBFR  cannot  possibly  reach  a  steady  state  in  which  there  is 
net  fusion  gain. 

In  several  of  their  progress  reports,  and  in  response  to  comments  by  ourselves  and 
by  A.  Carlson,11  RBM  argue  that  the  separation  between  the  azimuthal  rotation  velocities 
of  the  p  and  nB  beams  can  be  maintained  by  magnetic  forces,  arising  from  the  radial 
cross-field  drift  of  the  protons.  This  argument  is  not  viable.  If  the  proton  azimuthal 
velocity  is  assumed  to  be  at  steady  state,  with  the  proton-boron  friction  balanced  by  the 
vrxBz  magnetic  force,  then  it  is  necessary  that 


eVrB 

c 


=  VpBmpV<p  . 


(10) 


i.e. 


8.5  xlO-9 


nB  (cm  3) 
B  (gauss) 


cm/sec , 


(11) 


where  B  is  the  magnetic  induction  and  Qp  is  the  proton  cyclotron  frequency.  Taking  B  = 
105  G,  a  characteristic  value  for  the  RBM  scenarios,  we  find  that  the  required  radial 
velocity  is  Vr  =  8.5xlO~14nB  cm/sec.  Over  a  fusion  time  TF=5xl014nB~1  sec,  a  proton 
radial  velocity  Vr  of  this  magnitude  would  lead  to  radial  proton  motion  43  cm.  Since  the 
proton  ring  is  assumed  by  RBM  to  be  only  a  few  cm  wide  in  all  models  of  the  CBFR,  this 
would  constitute  radial  loss  of  protons  at  a  rate  considerably  exceeding  the  fusion  rate. 

O 

This  loss  mechanism  was  also  mentioned  by  Goldston  in  his  letter  to  Science. 


F.  Effect  of  Boron  Temperature 

Up  to  this  point,  we  have  assumed  that  the  boron  beam  is  cold.  But  collisional 
interactions  will  heat  the  boron  beam  as  well.  Scattering  of  protons  off  warm  nB  differs 
only  in  that  a  proton  energy  diffusion  term  is  also  present  in  the  Fokker-Planck  equation, 
which  in  effect  increases  the  thickness  of  the  spherical  shell  on  which  the  proton  velocity 
distribution  lies.  This  is  a  mechanism  of  temperature  equilibration  between  species, 
operating  in  addition  to  the  proton  heating  effects  previously  discussed,  which  are  due  to 
the  streaming  energy  between  the  beams.  Thus  effects  due  to  finite  boron  temperature 
only  add  to  the  proton  temperature,  and  further  detune  the  resonance. 
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G.  Proton-Proton  Scattering 


As  we  have  seen,  the  effect  of  p-nB  scattering  on  the  protons  is  primarily  to 
isotropize  the  protons  about  the  boron  velocity  Vb,  and  somewhat  more  slowly  to  extract 
energy  from  all  protons  at  the  same  rate.  Thus  on  the  rapid  time  scale  tpb,  p-nB 
collisions  spread  the  protons  out  on  an  energy  shell  with  characteristic  velocity  spread  vp 
~  IVp-VbI,  but  leave  all  of  the  protons  with  pretty  much  the  same  energy,  with  respect  to 
the  boron.  However,  we  must  also  consider  proton-proton  scattering,  which  drives  the 
proton  distribution  toward  an  isotropic  Maxwellian,  centered  on  the  mean  proton  velocity 
Yp,  and  with  thermal  velocity  vp.  The  proton-proton  momentum  scattering  rate  is 
vpp=7.2xl0‘13np(10  keV/Tp)372  ;  it  may  be  larger  or  smaller  than  the  p-nB  collision  rate 
depending  on  Tp  and  (np/nB),  but  vpp  is  orders  of  magnitude  larger  than  the  fusion  rate  Vf 
in  all  CBFR  designs  considered  by  RBM.  If  vpp  <  vpb,  proton-proton  scattering  modifies 
the  spherical-shell  proton  velocity  distribution,  by  filling  in  and  spreading  out  the 
constant-energy  shell,  and  over  a  longer  time  scale  vpp_1  driving  it  toward  Maxwellian.  If 
vpp »  vPB,  which  is  the  more  usual  case  for  CBFR,  proton-proton  collisions  keep  the 
proton  distribution  close  to  a  drifted  Maxwellian  at  all  times,  with  proton-boron  collisions 
steadily  increasing  Tp.  Since  both  vpp  »  vF  and  vpB  »  vF,  the  proton  streaming  energy 
will  in  all  cases  thermalize  long  before  the  fusion  time  vF_1,  thereby  totally  obliterating 
the  energy  resonance. 

H.  Ion-Electron  Collisions 

In  the  CBFR  scenario,  the  proton  kinetic  energy  in  the  electron  rest  frame,  Epe= 
(nip/2)(Vp-Ve)2,  is  much  larger  than  the  electron  temperature  Te,  but  the  electron  thermal 
velocity  ve  is  much  larger  than  the  proton-electron  relative  streaming  velocity  IVp-Vel.  In 
this  regime,  the  effect  of  ion-electron  collisions  on  the  ions  is  very  simple:  the  Fokker- 
Planck  equation  for  the  protons  reduces  to  simply  a  friction 

F=Vpe(v— Ve)  (12) 

on  the  protons.  The  diffusion  terms  in  the  Fokker-Planck  equation  are  smaller  by  a  factor 
Te/Epe,  and  can  be  neglected.  The  collision  frequency  for  this  process  is  Vpe=2.6xl0-14(10 
keV/Te)  ne,  independent  of  the  velocity  of  a  particular  ion.  These  well-known  facts  are 
proven  explicitly  in  the  Appendix.  Thus  the  effect  of  electron  scattering  on  the  protons  is 
to  draw  each  proton’s  velocity  toward  the  mean  electron  velocity  Ve.  Similarly,  the  effect 
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of  electron  scattering  on  the  boron  is  to  draw  each  11 B  nucleus’s  velocity  toward  the  mean 
electron  velocity  Ve,  in  this  case  at  a  rate  Vee  which  is  larger  than  Vpe  by  a  factor 
(Zb2/Ab)=25/1  1.  These  rates  are  faster  than  the  fusion  rate  Vp  for  most  CBFR  scenarios, 
and  for  those  cases  where  Te  is  not  very  large,  or  where  ne »  nB,  they  can  be  orders  of 
magnitude  larger. 

Let  us  consider  the  effects  of  these  collisions.  The  electron  inertia  is  always  small 
compared  to  the  ion  inertia,  and  thus  if  there  were  only  one  ion  beam  (say  the  proton 
beam),  the  electrons  would  very  quickly  be  speeded  up  to  the  proton  velocity  Vp,  with 
very  little  influence  on  the  protons.  However  the  electrons  are  also  subject  to  collisions 
with  the  nB  and  thus  they  settle  down  to  a  fluid  velocity  Ve  intermediate  between  Vp  and 
Vb.  Let  us  assume  that  Ve  is  in  a  steady  state.  The  effect  of  p-e  collisions  on  the  protons 
is  thus  to  draw  each  proton  toward  velocity  Ve.  Similarly,  nB-e  collisions  draw  each  UB 
nucleus  toward  Ve.  Thus  the  electrons  act  as  an  intermediary  that  couples  the  protons  and 
boron  through  a  friction  characterized  by  the  collision  frequencies  Vpe  and  Vse-  This 
friction  is  in  addition  to  the  friction  due  to  direct  p-nB  collisions.  In  most  of  the  CBFR 
scenarios  that  have  been  discussed  by  RBM,  the  electron-mediated  friction  is  fast 
compared  to  the  fusion  rate,  and  thus  effectively  prevents  maintenance  of  the  required 
beam  equilibrium.  In  the  scenario  presented  in  the  Science  article1,  where  np  =  400nB  and 
thus  ne  =  405nB,  and  Te  =  20  keV,  we  find  the  rates  to  be  Vpe  =  9xl0_15ne  =  1.8x103Vf  and 
Vse  =  2.  Ixl0-14rie  =  4.2x103Vf;  the  beams  are  virtually  stopped  in  their  tracks  by  electron 
friction. 

RBM  have  argued1  that  proton-electron  collisions  lead  to  significant  reduction  of 
the  proton  beam  temperature.  It  is  true  that  electron  collisions  can  have  a  cooling  effect 
on  the  protons,  but  this  is  merely  an  incidental  effect  of  the  friction  that  the  electrons 
exert  on  the  protons:  every  proton  is  drawn  toward  Ve,  and  thus  according  to  Eq.  (12)  the 
spread  in  the  proton  velocity  distribution  (as  seen  from  the  electron  mean  velocity  Ve) 
contracts  by  the  same  factor  as  the  reduction  in  IVp-Vel.  If  the  CBFR  is  to  succeed, 
parameters  must  be  chosen  so  that  Vpe  <  vF ;  otherwise,  proton-electron  friction  will  stop 
the  proton  beam.  However,  heating  of  protons  due  to  p-nB  pitch-angle  scattering 
proceeds  at  the  frequency  vpb,  which  we  have  seen  to  be  much  faster  than  Vp.  Thus, 
electron  cooling  cannot  prevail.  Additionally,  it  should  be  noted  that  for  very  non- 
thermal  velocity  distributions,  “cooling”  in  the  sense  of  reducing  the  spread  in  v-Ve  is  not 
necessarily  equivalent  to  reduction  of  the  spread  in  the  relative  proton-boron  velocity 
v-Vb,  nor  to  reduction  of  the  temperature  Tp  defined  in  Eq.  (2).  Consider,  for  example,  a 
proton  distribution  spread  on  a  constant-lv-Vsl  shell,  such  as  would  be  produced  by  p-B 
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collisions  alone,  as  seen  in  Fig.  3b.  Since  Ve  *  Vb,  the  electron  friction  draws  each  point 
on  the  proton  shell  toward  a  different  velocity  than  the  boron  friction.  As  a  result,  the 
shell  is  distorted  so  that  it  is  no  longer  a  spherical  shell  centered  on  Vb.  This  effect 
increases  the  spread  in  proton  energies  as  seen  by  the  boron  beam,  i.e.  in 


As  a  final  comment  on  the  subject,  we  note  that  in  most  cases  it  is  necessary  that 
the  electrons  be  somehow  maintained  at  a  large  temperature  Te,  in  order  to  reduce  the 
friction  rate  Vpe.  For  example,  in  their  most  recent  progress  report,14  RBM  propose  a 
scenario  in  which  Te  =  165  keV.  When  protons  scatter  off  warm  electrons,  a  proton 
energy  diffusion  term  is  present  in  the  Fokker-Planck  equation  in  addition  to  the  friction 
term,  which  increases  Tp  in  addition  to  the  effects  previously  mentioned. 

Collisional  Effects:  Summary  Statement 

In  Sections  3A-H  we  have  listed  at  least  eight  separate  dissipative  effects,  due  to 
ordinary  Coulomb  collisions,  each  of  which  detunes  the  resonance  and/or  destroys  the 
“colliding  beam”  equilibrium,  on  a  time  scale  fast  compared  to  the  fusion  time  scale. 
These  processes  are: 

1.  Isotropization  of  the  proton  beam  due  to  p-nB  scattering. 

2.  Loss  of  protons  scattered  into  v<p  <  0. 

3.  “Ohmic”  heating  of  the  protons,  a  consequence  of  any  force  that  maintains  the 
required  p-nB  streaming  velocity  in  the  presence  of  p-nB  scattering. 

4.  Radial  escape  of  protons,  if  the  radial  drift  velocity  is  assumed  large  enough  to 
maintain  the  required  p-nB  streaming  velocity. 

5.  Detuning  of  the  resonance  due  to  boron  acceleration,  resulting  from  nB-p  friction. 

6.  Detuning  of  the  resonance  due  to  proton  energy  loss,  resulting  from  p-nB  friction. 

7.  Stopping  of  the  relative  streaming  of  the  p  and  nB  beams,  due  to  electron  friction. 

8.  Thermalization  of  the  proton  momentum  spread,  due  to  p-p  collisions. 

Each  of  these  processes  increases  the  entropy.  For  successful  operation  of  the  p- 
nB  CBFR,  additional  collisional  or  other  effects  would  have  to  be  found  which  cancel,  to 
within  better  than  90%  (and  in  some  cases  99%),  all  of  these  effects  of  the  momentum 
and  energy  exchange  collisions.  Eight  mutually  consistent  miracles  is  a  lot  to  ask  for, 
even  in  the  fusion  business! 
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In  their  letters  to  Science ,  Nevins,10  Carlson,11  and  Rider  12  all  point  out  that  the 
protons  and  nB  cannot  be  maintained  as  cold  colliding  beams,  because  of  these  various 
ion-ion  scattering  processes.  In  their  response  to  Rider,24  Rostoker  and  Binderbauer 
appear  to  concede  this  point.  They  state: 

‘The  ion-ion  collision  term  is  certainly  a  problem,  particularly  if  the  investigation 
is  not  device  specific,  but  is  evaluated  with  generic  formulae.  There  are 
compensating  terms  for  the  CBFR. 

The  problem  of  ion-ion  collisions  can  be  avoided  with  a  different  mode  of 
operation  where  both  beams  have  the  same  average  velocity  and  temperature. 
Then  only  ion-electron  collisions  need  to  be  considered.” 

Obviously,  this  makes  no  sense  for  a  p-nB  reactor,  where  the  fusion  reaction  rate  is 
strongly  peaked  at  -600  keV,  the  energy  invested  in  accelerating  the  beams  exceeds  10% 
of  the  fusion  output  even  if  every  proton  fuses,  and  the  maximum  possible  thermal  fusion 
gain  is  only  a  little  larger  than  unity.  For  p-nB,  RBM  themselves  state6  that,  “it  is 
necessary  to  maintain  the  proton  and  boron  beams  at  an  average  energy 
difference  of  580  keV.  In  addition,  the  temperatures  of  the  beams  must  be 
substantially  less  than  1 40  keV.” 

4.  Heating  of  Protons  via  Proton-Boron  Collisions 

In  their  most  recent  ONR  progress  report,25  and  more  forcefully  in  a  widely 
circulated  letter  responding  to  our  criticism,26  RBM  concede  that  the  heating  or  cooling  of 
the  beams  is  a  crucial  issue,  and  then  argue  that  collisional  scattering  of  the  protons  by  the 
streaming  boron  beam  cools  the  protons,  i.e.  reduces  the  temperature  Tp  defined  by  Eq. 
(2).  In  their  response27  to  Nevins,10  RBM  claim  that  collisional  scattering  of  the  electrons 
by  the  streaming  boron  beam  cools  the  electrons.  These  are  amazing  claims.  It  should  be 
very  clear  that  scattering  a  cold  beam  of  light  particles  (protons  or  electrons)  off  a  cold 
beam  of  heavy  particles  (nB)  can  only  increase  the  temperature  of  the  light  beam.  The 
analogy  that  comes  to  mind  is  scattering  a  well-aligned  stream  of  tennis  balls  off  a 
distribution  of  bowling  balls.  Obviously,  in  the  frame  of  reference  in  which  the  bowling 
balls  are  stationary,  the  well-aligned  kinetic  energy  of  the  tennis  balls  is  quickly 
transformed  into  random  kinetic  energy,  i.e.  temperature. 

In  this  section  we  hope  to  lay  this  question  to  rest.  We  derive  an  explicit  equation 
for  the  time  evolution  of  the  proton  temperature,  based  on  the  Vlasov  /  Fokker-Planck 
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equation.  We  include  proton-boron  collisions,  proton-proton  collisions,  a  driving  electric 
field,  and  a  source  and  sink  of  protons.  The  calculation,  which  is  done  essentially 
exactly,  clearly  shows  that  the  effect  of  p-nB  collisions  is  to  heat  the  protons.  We  then 
retrace  the  derivation  offered  by  RBM,  to  elucidate  the  errors  that  led  them  to  conclude 
that  the  protons  are  cooled  by  p-nB  collisions. 

This  section  includes  more  mathematical  detail  than  most  of  the  paper.  Since 
proton  heating  has  been  a  point  of  controversy  between  us  and  RBM,  we  have  tried  to 
present  the  derivations  in  enough  detail  to  make  them  accessible  to  the  general  reader, 
and  hopefully  to  resolve  any  disagreements. 

Our  treatment  begins  with  the  kinetic  equation  controlling  the  normalized  proton 
velocity  distribution  fp(v,t), 


npeE  afp 
mp  dv 


'npav,Ffp+2lXp9v,av 


±vv± 


:D±fP+v^ 


pp 


(13) 


+ hj5(v-  V! )  -  vpnp8(v- v2) . 


The  second  term  on  the  LHS  is  the  standard  Vlasov  term  representing  the  effect  of  an 
electric  field.  We  have  in  mind  that  E  is  primarily  an  azimuthal  electric  field  E<p  which 
may  be  present  to  drive  the  proton  streaming.  The  first  two  terms  on  the  RHS  are  the 
Fokker-Planck  friction  and  diffusion  terms  representing  p-nB  scattering,  as  discussed  in 
the  Appendix.  These  terms  are  the  focus  of  our  interest  here.  We  note  that  the  diffusion 
tensor  D  is  purely  transverse  for  scattering  off  a  cold  boron  beam,  as  proven  in  the 
Appendix.  We  have  also  included,  for  completeness,  a  collision  integral  representing  p-p 
collisions,  although  this  term  contributes  nothing  to  the  fluid  equations.  The  last  two 
terms  represent  the  source  of  new  protons  (the  injected  proton  beam),  and  the  sink  of 
protons  due  to  fusion  events.  These  terms  are  also  included  for  completeness,  although 
they  make  only  a  small  contribution  to  the  proton  temperature  equation.  The  quantity 
vF  is  an  average  fusion  rate,  which  will  necessarily  be  somewhat  smaller  than  the  peak 
fusion  rate  Vf  (at  energy  Eo)  defined  previously.  We  assume  that  the  protons  are  injected 
at  a  well-defined  azimuthal  velocity  vi  (a  very  realistic  assumption),  and  that  fusion 
reactions  occur  at  a  well-defined  velocity  \2  (an  unrealistic  assumption  that  is  made  only 
to  simplify  the  presentation).  Our  treatment  assumes  that  both  beams  are  uniform  in  the 
azimuthal  direction  cp  (the  streaming  direction).  We  also  neglect  inhomogeneities  in  the 
radial  and  axial  directions,  as  they  are  not  relevant  to  the  present  discussion.  A  more 
rigorous  formulation  would  express  the  results  as  a  radial  average,  as  is  done  by  RBM  in 
their  analysis.  This  complicates  the  notation,  but  does  not  change  the  results. 
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We  proceed  to  derive  fluid  equations  by  taking  velocity  moments  of  Eq.  (13). 
This  procedure  can  be  carried  out  rigorously,  because  the  beams  have  been  taken  to  be 
homogeneous.  We  first  take  the  velocity  integral  of  Eq.  (13),  to  arrive  at  a  density 
continuity  equation.  This  equation  is  trivial: 


dnp 

dt 


hi  -  vFn 


p  ‘ 


(14) 


In  steady  state,  the  beam  injection  rate  must  equal  the  fusion  loss  rate,  so  that  (14)  reduces 
simply  to 

dn„ 


dt 


=  0, 


ni  =  vFnp. 


(15) 


We  shall  assume  Eqs.  (15)  are  valid.  We  next  take  the  moment  Jd3v  mpv,  in  order  to 
derive  the  momentum  equation: 


m. 


^}d3vvfp+eE-J d3v-^-v  =  -mpJd3w^  Ffp  +  vFmp(Vj  -v2).  (16) 


The  Fokker-Planck  diffusion  term  and  the  p-p  collision  term  make  no  contribution,  as 
these  terms  conserve  total  proton  momentum.  Using  integration  by  parts,  Eq.  (16) 
reduces  exactly  to 

dV_ 

mp  — =  eE  +  mpjd3vFfp+vFmp(v,-v2)  .  (17) 

Finally,  we  take  the  moment  Jd3vymp(v-  Vp)2  to  derive  a  temperature  equation. 

The  p-p  collision  term  makes  no  contribution,  since  this  term  conserves  total  proton 
energy.  Using  the  definition  (2)  of  Tp,  and  performing  integration  by  parts  for  the  Fokker- 
Planck  terms  (but  no  approximations),  the  remaining  terms  can  be  put  in  the  form 

2^=eEJdM*-vph +mpldM*-vp)  FfP +impJ‘*svDfP 

(18) 

+  vFmp(v1-v2)^^iY^-Vpj, 


where  D  is  the  trace  of  the  Fokker-Planck  diffusion  tensor.  We  note  that  the  first  term  on 
the  RHS  of  (18)  integrates  to  zero,  i.e.  the  electric  field  makes  no  direct  contribution  to 
dTp  /  dt.  We  now  use  the  explicit  forms  for  F  and  D_l  obtained  in  Eqs.  (A10)  of  the 
Appendix: 
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F  =  -vpb(v--V b). 


(19a) 


D  =  2v 


nir 


pB 


mB+mp 


(v-VB)2, 


(19b) 


where  vpB  is  the  proton-boron  collision  frequency, 
4rcnBe4ZB2  mB  +  mp  ^ 


m 


,(v  —  VB) 


(20) 


and  X  is  the  Coulomb  logarithm.  The  collision  frequency  vpB  is  a  slowly-varying  function 
of  |v-VBl,  but  lv-VBl  is  nearly  the  same  for  all  protons  in  the  cold  beam.  Thus  we  shall 
treat  vpB  as  a  constant.  (This  very  reasonable  approximation  is  also  used  by  RBM  in  all 
of  their  calculations.)  The  momentum  equation  (17)  then  reduces  exactly  to 

mp  ^  =  eE  -  mpvpB(vp  -  VB)  +  vFmp(v,  -  v2 ) ,  (21) 


while  the  temperature  equation  (18)  becomes 


3dTr 
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Collecting  terms,  Eq.  (22)  reduces  to 
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— V„ 


.  (23) 


The  first  term  of  Eq.  (23)  represents  heating  of  the  protons  due  to  p-nB  scattering.  The 
CBFR  requires  mp(Vp- VB)2  ~  1.2  MeV  »  Tp.  Thus  the  second  term  on  the  RHS  is  a 
small  correction,  and  it  is  clear  that  the  scattering  process  results  in  proton  heating,  not 
cooling.  The  last  term  is  also  a  small  correction,  since  the  fusion  rate  vF  is  much  smaller 
than  vpB  and  V1-V2  may  be  expected  to  be  less  than  IVp-VBl. 

We  now  retrace  through  Ref.  14  to  try  to  understand  how  RBM  arrive  at  the 
contrary  conclusion  that  proton-boron  scattering  can  cool  the  protons.  The  derivation 
begins  in  Appendix  A  of  Ref.  14,  on  p.  9,  and  then  continues  in  the  main  body  of  the 
paper,  on  p.  4.  In  Appendix  A,  the  discussion  is  in  terms  of  the  effects  of  scattering  on 
electron  temperature,  but  on  p.  4  the  previous  equations  are  applied  to  proton 
temperature.  The  justification  for  some  of  the  steps  in  the  derivation  is  not  evident. 
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The  RBM  derivation  begins  with  Eq.  (4),  p.  9,  in  Appendix  A  of  Ref.  14,  which 
we  restate  in  our  notation  as 


|^-  =  mpjd3v[(v-Vp)F  +  iD]fp(v), 


(24) 


[RBM  use  the  notation  (Av)  for  our  F,  and  (AvAv)  for  our  D.  We  have  rewritten  the 
equation  for  protons  instead  of  electrons,  and  we  have  suppressed  the  radial  averaging 
that  is  applied  to  every  term  and  makes  the  notation  look  much  more  complicated.] 
Equation  (24)  is  similar  to  our  Eq.  (18).  The  source  and  sink  terms  are  not  included,  but 
they  are  not  relevant  to  the  present  discussion.  However,  RBM  then  continue  with: 


“We  note  that  if  np  and  Vp  are  constant, 


adTp 

2  dt 


•E  + 


mp|d3v[vF  +  ^D]fp(v).” 


(25) 


In  going  from  Eq.  (24)  to  Eq.  (25),  a  term  -Vp-F  has  been  replaced  by  eVp-E.  We  can 
only  speculate  that  the  intention  is  to  insert  the  effect  of  an  electric  field,  and  to  reason 
that  if  a  steady  state  exists  the  electric  field  must  take  the  value  eE  =  -F  [see  Eq.  (17)]  so 
as  to  counterbalance  the  friction  F.  Since  the  proton  flow  is  primarily  in  the  azimuthal 
direction,  E  and  F  would  also  have  to  be  vectors  directed  primarily  in  the  azimuthal 
direction.  If  this  substitution  is  made,  the  eVp  E  term  must  be  regarded  as  an  essential 
part  of  the  effect  of  p-nB  scattering  on  dTp/dt.  However,  this  term  is  simply  ignored 
throughout  the  remainder  of  RBM’s  Appendix  A,  up  to  the  end  of  p.  10.  The  discussion 
is  concerned  entirely  with  the  evaluation  of  the  second  term  on  the  RHS  of  Eq.  (25), 
which  RBM  seem  to  regard  as  a  complete  description  of  the  p-nB  interaction.  In  this 
way,  RBM  finally  arrive  their  Eq.  (12),  p.  10,  which  is  (in  our  notation) 

i^  =  -AvpBmpVB(Vp-VB)  (26) 

where  A  is  a  coefficient  that  includes  area-weighting.  This  is  essentially  the  equation  that 
is  subsequently  quoted  in  the  main  text  of  their  paper  [Eq.  (9),  p.4],  and  in  subsequent 
letters,  as  indicating  that  the  effect  of  p-nB  collisions  is  to  cool  the  protons. 

We  do  not  understand  the  reasoning  behind  these  steps,  which  do  not  constitute  a 
correct  treatment  of  the  steady  state  nor  of  the  electric  field  effect,  and  which  are  clearly 
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inconsistent  with  the  derivation  of  Eqs.  (18,22,23)  directly  from  the  Vlasov/Fokker- 
Planck  equation  (13).  It  should  be  clear  just  from  its  form  that  Eq.  (26)  cannot  be  correct: 
Temperature,  as  defined  in  Eq.  (2)  -  the  same  definition  used  by  RBM  -  is  a  quantity  that 
is  independent  of  the  frame  of  reference  of  the  calculation.  Equation  (26)  does  not  have 
this  property:  on  the  RHS,  the  factor  (Vp  -  VB)  is  frame-invariant,  but  the  factor  VB  is 
not.  By  merely  doing  the  calculation  in  a  different  (rotating)  frame  of  reference,  the 
“cooling”  effect  can  be  converted  to  a  heating  effect,  of  arbitrary  magnitude.  If  RBM  had 
followed  through  with  the  calculation  directly  from  Eq.  (24)  or  (25),  they  would  have 
arrived  at  an  equation  essentially  equivalent  to  our  Eq.  (23).  Indeed,  they  themselves, 
immediately  after  writing  Eq.  (26),  comment  (at  the  top  of  p.  11,  Ref.  14)  that 

“If  Ve  was  not  assumed  to  be  constant  the  result  would  be 

^  =  -AvpBmp(Vp-VB)2.”  (27) 

Equation  (27),  which  is  equivalent  to  our  Eq.  (23),  does  have  the  correct  property  of 
frame-invariance,  and  is  the  correct  result  in  steady  state  or  in  a  transient  situation. 

Strangely  enough,  RBM  do  eventually  return  to  the  term  eVp-E  in  Eq.  (25),  but 
only  to  include  the  radial  component  eVprEr,  which  is  associated  with  the  slow  radial 
diffusive  drift  of  the  protons  in  a  steady  state  burning  plasma.  This  contribution  is 
discussed  on  p.  11  of  RBM’s  Appendix  A,  and  also  on  pp.  4-5  of  their  main  text.  As 
noted  by  RBM,  this  radial  term  is  small. 

5.  Other  Issues  Related  to  the  p-11B  CBFR 

In  this  section,  we  comment  very  briefly  on  some  additional  technological  and 
scientific  issues.  Although  these  issues  do  not  have  the  force  of  obvious  show-stoppers, 
they  do  appear  to  raise  serious  concerns  regarding  the  viability  and  utility  of  a  CBFR 
reactor,  particularly  in  regard  to  the  naval  shipboard  power  application,  where 
compactness  is  of  the  essence. 

A.  Beam  Technology  and  Size  Constraints 

The  ONR  reports3  of  RBM  are  largely  concerned  with  the  application  of  this 
fusion  scheme  as  a  power  source  for  naval  ships.  For  this  application,  minimization  of 
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device  size  is  crucial,  and  RBM  devote  considerable  effort  to  this  issue.  In  the  Science 
article  as  well,1  continued  interest  in  the  naval  application  is  evident,  and  indeed  Fig.  3  of 
Ref.  1  is  intended  to  show  the  compact  nature  of  a  naval  reactor.  However,  nowhere  is 
there  any  discussion  of  the  proton  beam  source,  which  is  simply  sketched  in  that  figure  as 
a  box  approximately  Ix2mx2min  size. 

For  purely  steady  state  operation,  neutral  beam  injection  would  be  required.  For 
the  naval  application,  the  RBM  scheme  requires  20  amps  of  neutral  H  at  energy  >  600 
keV,  i.e.  an  input  beam  power  of  >  12  MW,  to  produce  a  maximum  fusion  power  of  about 
160  MW  if  every  proton  fuses.  The  preeminent  existing  program  in  neutral  H  beam 
development  is  being  conducted  in  Japan,  for  the  purpose  of  heating  the  JT  60-U 
Tokamak.28  This  facility  currently  delivers  2.5  MW  of  H  atoms  at  350  keV,  and  planned 
upgrades  will  bring  it  to  10  MW  at  500  keV,  comparable  to  the  beam  required  for  the 
CBFR.  However,  the  volume  required  for  neutral  beam  production  and  focusing  is 
enormous,  compared  to  the  scale  size  of  2m  sketched  by  RBM  in  Fig.  3  of  Ref.  1.  The 
sources  are  approximately  6m  long,  and  because  the  beams  are  generated  with  very  large 
diameter  and  must  be  focused  down  with  good  collinearity  (5  mrad),  the  focusing  length 
is  approximately  18  m.  There  is  no  known  neutral  beam  source  technology  which  is  close 
to  the  compactness  needed  for  the  naval  application. 

A  more  fundamental  concern  is  that  trapping  of  the  neutral  H  beam  in  the  CBFR 
plasma  would  appear  to  be  very  inefficient,  because  the  proton  ring  envisioned  by  RBM 
presents  a  very  small  target.  To  be  trapped  in  the  plasma,  a  beam  H  atom  must  be  ionized 
by  interaction  with  the  plasma.  In  the  regime  of  interest,  electron-impact  is  the 
predominant  ionization  mechanism,  with  an  ionization  cross-section  of  about  10  cm 
at  Te  =  20  keV.  (The  charge-exchange  cross-section  is  even  smaller  at  proton  energy  600 
keV;  extrapolation  from  Fig.  3.1  of  Brown30  indicates  that  it  is  about  5xl0’20  cm2.)  For 
electron  density  ne  =  2xl015  cm-3,  the  mean  free  path  for  ionization  of  an  H  atom  is  thus 
~500  cm.  But  for  injection  perpendicular  to  the  magnetic  field,  the  CBFR  originally 
proposed  with  radius  r  =  30  cm  and  thickness  Ar  =  4  cm  presents  a  target  plasma  at  most 
30  cm  long  along  any  chord,  as  shown  in  Fig.  4,  and  15  cm  on  the  average.  Some  of  the 
more  recent  designs  by  RBM3  envision  a  proton  ring  with  radius  up  to  80  cm.  Even  if  Ar 
=  10  cm,  the  average  chord  is  only  40  cm,  so  that  less  than  10%  of  the  beam  will  be 
trapped.  Furthermore,  the  narrowest  beam  that  can  be  produced  by  the  JT  60-U  facility 
(or  any  planned  follow-on  technology)  is  much  wider  than  the  maximum  acceptance  of 
the  CBFR  plasma.  Thus  neutral  beam  injection,  using  any  reasonable  extrapolation  of 
today’s  technology,  would  not  appear  to  be  a  viable  driver  for  the  p-nB  CBFR. 
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An  alternative  approach  is  repetitively  pulsed  injection  of  a  neutralized  ion  beam 
(i.e.  an  ion  beam  accompanied  by  an  equal  and  opposite  charge  density  of  electrons),  with 
parameters  chosen  so  that  the  ion  plasma  frequency  exceeds  the  ion  gyrofrequency.  Then 
the  neutralized  ion  beam  can  cross  the  vacuum  magnetic  field  via  the  ExB  drift,  where  E 
is  the  polarization  field.  The  polarization  field  is  shorted  out  in  the  dense  plasma,  and  the 
beam  is  trapped.  The  most  powerful  present-day  repetitively  pulsed  ion  beam  sources, 
e.g.  the  RHEPP-II  ion  accelerator  at  Sandia,32  produce  a  time-averaged  current  of  only  0.1 
amp,  two  orders  of  magnitude  below  the  CBFR  requirement  for  a  100  MW  naval  reactor, 
and  a  time-averaged  power  of  300  kW,  a  factor  of  40  below  the  CBFR  requirement. 
(However,  the  design  voltage  of  2.5  MV  is  three  to  four  times  greater  than  the 
requirement  for  CBFR.)  This  source  occupies  a  volume  of  about  400  m3,  over  two  orders 
of  magnitude  larger  than  the  ion  accelerator  sketched  in  Fig.  3  of  Ref.  1.  The  most 
advanced  proposed  repetitively  pulsed  electron  sources,  e.g.  the  E-Scrub  electron 
accelerator33  proposed  by  SAIC  (for  which  minimization  of  size  was  a  prime  objective), 
would  occupy  a  volume  of  500  m3  for  12  MW  of  average  power  at  voltage  800  kV. 
Electron  sources  are  typically  considerably  more  compact  than  ion  sources.  There  is  no 
known  technology  which  can  supply  the  rep-rated  ion  beams  for  the  CBFR  with  the 
compactness  envisioned  for  the  naval  application. 

B.  Experimental  Evidence  on  Field-Reversed  Ion  Rings 

There  is  a  lengthy  publication  record  on  field-reversed  configurations.34’35 
However,  this  body  of  work  is  concerned  primarily  with  plasma-filled  RFC’s.  The 
experimental  evidence  on  RFC’s  powered  by  injected  ion  rings  is  scant,  and  it  is  not 
encouraging  in  regard  to  the  prospects  for  maintaining  a  thin  ring  at  the  resonant  energy 
necessary  for  the  p-nB  reaction.  Experiments  were  done  in  the  1980’s  at  NRL36  and 
Cornell37.  In  both  cases,  the  ion  ring  spread  out  radially  so  that  the  thickness  Ar  became 
comparable  to  r.  These  experiments  were  not  specifically  designed  to  produce  thin  rings, 
but  they  do  suggest  that  maintaining  a  thin  ring  may  be  very  difficult. 
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C.  Instabilities 


As  RBM  noted  in  their  response  to  Rider,24  “The  list  of  possible  instabilities  is 
endless.”  However,  we  shall  simply  comment  here  that,  on  general  theoretical  grounds 
one  might  expect  that  a  thin  field-reversed  ion  ring  equilibrium,  if  it  could  somehow  be 
produced,  would  be  susceptible  to  major  disruptions.  After  nearly  40  years  of  tokamak 
research,  the  problem  of  major  disruptions  still  has  not  been  solved.  Tokamaks  are  low- 
beta  devices,  so  it  is  only  the  poloidal  field  energy  which  is  released  in  a  disruption.  In 
the  CBFR,  virtually  all  of  the  field  energy  could  be  released.  For  the  30  cm  radius  beam 
with  a  40  cm  inner  wall,  the  magnetic  free  energy  of  the  field  reversed  configuration  (i.e., 
the  difference  between  the  magnetic  energy  in  the  properly  configured  CBFR,  and  the 
magnetic  energy  of  a  uniform  solenoidal  field  with  the  same  flux)  is  about  20  MJ 
(equivalent  to  about  5  pounds  of  TNT)  per  meter.  This  far  exceeds  the  energy  content 
(about  4  MJ  per  meter)  of  the  thin  ion  ring  which  separates  the  reversed  field  regions.  A 
single  major  disruption  on  this  scale  would  most  likely  destroy  the  superconducting 
magnet. 

6.  Conclusions 


The  use  of  the  p-nB  fusion  reaction  is  alluring,  because  of  the  absence  of  neutrons 
and  the  possibility  of  direct  conversion  of  fusion  energy  to  electricity.  These  would  be 
major  advantages  in  a  civilian  power  reactor,  and  they  would  be  overwhelming 
advantages  in  a  naval  power  source.  Unfortunately,  as  has  been  noted  by  many 
authors,5’8’9  the  energy  balance  for  p-nB  fusion  is  so  marginal  that  virtually  everything 
must  work  out  perfectly  to  make  net  gain  possible.  In  particular,  to  achieve  substantial 
gain  it  would  appear  to  be  necessary  to  operate  with  colliding  p  and  nB  beams  (as 
proposed  by  Rostoker,  Binderbauer  and  Monkhorst)  rather  than  with  a  thermal  plasma,  in 
order  to  exploit  the  resonance  in  the  fusion  cross-section  at  about  600  keV. 

Unfortunately,  for  p-nB  there  are  many  Coulomb  collisional  processes  which  are 
orders  of  magnitude  faster  than  the  fusion  rate.  We  have  pointed  out  eight  processes, 
each  of  which  has  the  effect  of  dissipating  the  cold-colliding-beam  equilibrium  by  heating 
the  beams  and/or  bringing  their  velocities  together.  There  is  no  possibility  of  preserving 
the  required  equilibrium  for  long  enough  to  allow  a  substantial  fraction  of  the  injected 
protons  to  undergo  fusion,  and  therefore  no  possibility  of  achieving  net  gain  from  the 
CBFR.  These  arguments  would  appear  to  apply  as  well  to  any  type  of  colliding-beam 
scheme  for  p-nB  fusion. 
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We  have  additionally  pointed  out  that  there  is  no  reason  to  expect  that  the  CBFR 
concept  could  be  implemented  in  a  compact  form  suitable  for  a  naval  power  source, 
inasmuch  as  the  technology  for  producing  and  injecting  beams  involves  very  large 
equipment.  Some  other  difficulties,  related  to  beam  trapping  efficiency  and  to  plasma 
stability,  have  also  been  discussed  in  passing. 

This  work  was  supported  by  ONR. 
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Appendix:  Analysis  of  the  Fokker-Planck  Equation 


The  Rosenbluth-MacDonald-Judd  form  of  the  Fokker-Planck  equation,38  which 
represents  the  effect  of  multiple  Coulomb  collisions  in  a  plasma,  is 
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where  F(v)  and  D(v)  are  the  dynamical  friction  and  velocity  diffusion  tensor,  given  by 
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j  is  the  species  whose  distribution  function  fj(v)  is  being  followed,  k  is  the  scatterer 
species  [Eqs.  (A2)  should  be  summed  over  various  scatterer  species],  and  X  is  the 
Coulomb  logarithm.  If  the  scatterer  distribution  is  isotropic  about  its  own  mean  velocity 
yk,  we  have  shown39  that  F  and  D  can  be  reduced  to  the  following  forms: 
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where  u=v-Vk,  Du  is  the  parallel  diffusion  (the  diagonal  component  of  D  along  u),  and 
Di  is  the  transverse  diffusion  (the  diagonal  component  of  D  perpendicular  to  u). 
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Scattering  of  Energetic  Protons  by  Electrons 


For  proton  scattering  off  electrons,  the  friction  coefficient  F  is  multiplied  by  a 
very  large  factor  mp/me,  which  is  lacking  in  the  diffusion  coefficients.  This  leads  one  to 
expect  friction  to  be  dominant  over  diffusion.  Looking  at  the  situation  more  carefully,  it 
becomes  clear  that  friction  dominates  over  diffusion  if  the  proton  kinetic  energy  is  large 
compared  to  the  electron  temperature.  The  proof  of  this  is  as  follows.  With  velocities 
referenced  to  the  electron  fluid  velocity  Ve,  the  fractional  change  in  proton  velocity  in  a 
time  interval  At  is 
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whereas  the  fractional  spread  in  proton  kinetic  energy  due  to  diffusion  is 
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Evaluating  Eqs.  (A3),  we  find  that  in  the  limit  of  very  high  proton  energy,  where  the 
proton  velocity  u  exceeds  the  electron  thermal  velocity  ve, 


Au  /  u  _  mP 
Au2  /  u2  2me 


(A5) 


and  in  the  intermediate  regime  where 
2  mpu  »  Te ,  but  u  <  ve, 

we  find 

Au  /  u  mu2 

T  2~.  2  =  ~ or  ^  1  • 

Au  /  u  2Te 


(A6) 


(A7) 


Thus,  for  scattering  of  energetic  protons  (i.e.  protons  with  mu2  »  2Te)  by  electrons,  the 
diffusion  terms  in  the  Fokker-Planck  equation  are  always  small,  and  the  Fokker-Planck 
equation  reduces  to  simply  a  frictional  drag  on  the  protons.  Furthermore,  in  the 
intermediate  proton  velocity  regime  where  (A6)  holds,  (A3 a)  can  integrated 
approximately,  to  give 
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We  see  that  in  this  regime  the  friction  is  proportional  to  proton  velocity  u,  i.e.  the 
collision  frequency  Vpe  is  velocity-independent. 


Scattering  by  a  Cold  Scatterer  Species 


If  the  scatterer  species  k  is  cold,  i.e.  has  distribution  function 


fk(u)  =  S(u), 

k  4rcv2 

Eqs.  (A3)  can  be  evaluated  trivially.  The  result  is 
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For  cold  scatterers,  the  parallel  diffusion  vanishes,  and  the  Fokker-Planck 
equation  reduces  to  friction  and  transverse  diffusion.  In  the  limit  nij/mk  0,  there  is  no 
energy  exchange  during  an  elastic  collision,  and  therefore  one  expects  the  Fokker-Planck 
equation  for  fj(v)  to  reduce  to  diffusion  on  a  constant-energy  spherical  surface.  In  the 
more  general  case  where  the  scatterers  are  cold,  but  it  is  not  assumed  that  mj/mic  is  small, 
one  expects  the  Fokker-Planck  equation  to  reduce  to  diffusion  on  a  constant-energy 
spherical  surface,  plus  self-similar  contraction  of  the  spherical  surface  due  to  friction.  We 
shall  now  prove  that  this  is  correct,  by  rewriting  Eqs.  (Al,10)  in  spherical  coordinates 
(u,0,\|/),  centered  on  the  particle  velocity  u,  with  0  the  polar  angle  and  the  azimuthal 
angle.  Noting  that  Dj.  is  the  only  nonvanishing  component  of  the  tensor  D,  the  second 
term  on  the  RHS  of  (Al)  can  be  written 
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Since  Dj.  from  Eq.  (A  10c)  is  independent  of  0  and  (p,  this  reduces  to 
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According  to  Eqs.  (A  10), 
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and  thus  we  can  combine  the  terms  of  Eq.  (Al)  to  get 
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The  first  term  of  (A14)  is  the  friction  term  in  spherical  coordinates.  Equation  (A15) 
shows  that  this  term  vanishes  in  the  limit  of  light  particles  scattering  off  heavy  particles. 
The  second  term  of  (A14)  is  just  diffusion  on  the  spherical  constant  energy  shell. 

Equation  (A14)  has  an  exact  solution  in  the  form 
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The  second  term  on  the  LHS  of  (A  18)  is  small  in  the  limit  of  light  particles  scattering  off 
heavy  particles.  If  the  distribution  fj(v)  initially  has  zero  energy  spread,  then  the  Fokker- 
Planck  equation  separates  into  diffusion  on  a  constant-energy  shell,  and  self-similar 
collapse  of  the  energy  shell  at  the  friction  rate  given  by  Eq.  (A18).  There  is  no  spread  in 
the  particle  energies  at  any  given  time.  This  conclusion  is  relevant  to  the  evolution  of  the 
proton  velocity  distribution  under  the  influence  of  p-nB  collisions,  since  both  the  protons 
and  the  boron  are  assumed  to  begin  as  cold  beams  with  very  little  energy  spread. 
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Figure  1.  Fusion  cross-section  for  unpolarized  p-nB,  as  a  function  of  collision  energy  E,  according  to 
Becker  et  al.17  (A  smaller  cross-section  is  given  by  Ref.  16.)  RBM  argue  that  the  Becker  cross-section  can  be 
further  increased  by  a  factor  of  1.6  for  100%  spin-polarized  fuel. 
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Figure  3.  (a)  Ion  velocity  distributions  for  the  cold-colliding-beam  equilibrium  assumed  for  the  CBFR.  The  dark  circle  at  the 
origin  represents  the  nB  velocity  distribution,  while  the  circle  to  the  right  represents  the  proton  velocity  distribution. 

(b)  Proton  velocity  distribution  resulting  from  p-nB  momentum  scattering,  at  time  ~v  g*1.  The  protons  have  diffused  over  the 
right  hemisphere  of  the  constant-u  sphere. 

(c)  Proton  velocity  distribution,  at  time  >vpB'\  if  p-nB  momentum  scattering  were  the  only  process  at  work.  The  proton 
diffusion  has  progressed  to  the  point  where  the  proton  distribution  is  uniform  on  the  constant-u  sphere. 
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Figure  4.  Geometry  for  neutral  H  atom  injection  into  the  CBFR.  To  be  trapped,  the  H  atoms  must  be  ionized 
while  traversing  the  chord  of  length  <Lmax. 
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